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Abstract. The determination of confidence intervals of stress orientation is a 
crucial element in the discussion of homogeneity or heterogeneity of the stress 
field under study. The error estimat es provided by the gr i d sea rch method 
Focal Mechanism Stress Inversion of iGephart and For syth ( 1984]) have been 
shown to be too wide but the reasons for this failure have escaped elucidation. 

Through the use of directional statistics and synthetic focal mechanisms, I 
show that the grid search methodology does yield appropriate uncertainty es- 
timates. The direct perturbation of the synthetic focal mechanisms introduces 
bias which leads to confidence intervals which become increasingly too wide as 
the amount of perturbation increases. The synthetic data also show at what 
point the method fails to overcome this bias and when confidence intervals will 
be too wide. The indirect perturbation of the focal mechanisms by perturbing 
the generating deviatoric stress tensor generates synthetic data devoid of bias. 
Inversion of these data sets yields correct confidence intervals. 

The Focal Mechanism Stress Inversion method is vindicated as a highly 
effective method, and with the use of appropriate directional statistics, its 
results can be assessed and homogeneity or heterogeneity of the stress field 
can be discussed with confidence. 



1. Introduction 

Gephart and Forsvtb ( 1984 ) proposed a grid search method of inverting focal 



mechanisms to obtain the stress tensor (focal mechanisms stress inversion, hence- 
forward FMSI) , in which stress field parameters are tried systematically against the 
focal mechanism orientations and the misfit calculated. They defined the misfit as 
the sum of the minimum angle needed to bring the slip direction of each focal mech- 
anism into line with the resolved shear stress on the fault plane. Both planes of each 
focal mechanism are tried and the smallest misfit serves to differentiate fault plane 
from auxiliary plane. They incorporated t he one-norm measur e as m isfit crite- 
rion, adopting the methodology proposed by IParker and McNutt (Il980l). Different 



techniques for such inversions had been proposed already, and I Gephart] (|1990bl ) 
discussed their relative merits, particularly in the context of focal mechanisms and 
th e problems asso ciated with the presence of two nodal planes 



Michaell (1987) proposed a different method, relying on a linearisation and boot- 



strapping to invert focal mechanisms for stress tensor calculation. He drew atten- 
tion to some differences in the size of the confidence intervals obtained by these two 
different methods. 
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It became quickly apparent that the discrepancies in confidence intervals had 
major repercussions on the study of actual stress fields. Different methods led 
different groups of researchers to different conclusions regarding the spatial and 
temporal variatio n in the stress regime in Southern California. This unsatisfactory 
state of affairs led Hardebeck and Hauksson ( 200ll) to test thoroughly the different 
methodologies using synthetic data. They showed that the confidence intervals 
estimated by FMSI were much too large, but they did not succeed in determining 
the reasons for such overestimates. 

Here, I attempt to unco ver the nature of the confi d ence i ntervals estimated by 
FMSI by investigating the lHardebeck and Hauksson! (|200ll ) analysis. Initial at- 
tempts to resolve the problem of these too wide confidence intervals brought to light 
some inac curacies in the derivation an d application of the one- norm measure as pro - 
pos ed bv IParker and McNutt 
and lHardebeck and Hauksson! 




and applied b y iGephart and Forsvtbl (|l984h 
(|Revetsl 12009). These corrections proved in- 



sufficient to lessen significantly the width of the confidence intervals (Hardebeck, 
pers. comm.). 

This left two avenues of investigation: the nature of the statistics used to calcu- 
late the confidence intervals, and the methodology of calculating the synthetic data 
used to test the method. 



2. Directional statistics 
Fisher ( 1953) pointed out that the initial development of the theory of errors 



by Gauss aimed to help astronomers and surveyors with their accurate angular 
measurements. Because of the accuracy of their measurements, the linear approxi- 
mation to which Gauss resorted was both appropriate and effective. However, when 
the errors become substantial, the linear approximation is no longer valid and the 
topological framework has to be taken into account. 

The spherical mean direction R and variance S are defined as 

N N N 

2 / \ -> \ 2 / \ -> \ 2 



(1) # = +(£"*) +(E 

i=l i=l 

and 

(2) S=(N-R)/N 



ii. 



wher e U, m%, ni are the direction cosines of the angular measurements (jMardia 



1972). 



The spheri cal equivalent of the normal distribution is the Fisher distribution 



(jFisherl 119531) . defined by 



(3) df = 



2 sinh k 



in which k is a measure of concentration. When k approaches zero, the distribution 
becomes uniform over the entire sphere whereas for k large, the distribution is 
confined to a small region of the sphere around the maximum. In the latter case, 
the distribution comes close to a two-dimensional (isotropic) normal distribution 
where k plays the role of the inverse of the variance. 
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The maximum likelihood estimate of k for the distribution on a sphere is for 
larger values of k 

il.5(«) 



(4) 



R 



I0.5 («) 



where I n is the modified Bcsscl function of the first kind and the n-th order ( Watson 



1944 ) and R the mean of R (jWatson and Williams! . Il956t iMardiai Il972h . For large 
values of n, the most likely value is given by 

N-l 

^ K =N^R 



(jFisherl . Il953t IWatsonl . Il960h . 



Confidence intervals for k can be determined from the relation 
(6) 2k(N-R) 
( Watson and Williamsl . Il95fit IMardiai Il972h . 



X2N-2 



3. Application to stress tensor inversion 

To calculate and test the confidence intervals of stress tensors inverted from focal 
mecha nisms, I adapted and modified the method used bv lHardebeck and Haukssonl 
(|200lh . taking particular care to adhere to the strictures and requirements of direc- 
tional statistics. 

3.1. Generating synthetic focal mechanisms. Sets of focal mechanisms are 
made up from spherically randomly oriented fault planes with slip directions de- 
termined by a chosen stress tensor. The azimuth of the fault planes is chosen 
from uniformly distributed random numbers in the [0, 2ir] interval, while the dip is 
randomly taken out of the [0, 1] interval through 

(7) di G arccos[0, 1] 

to avoid the (spherical) bias that the direct selection from the [0, n/2] interval would 
cause. Sets contain either 20, 50 or 100 fault planes. 

The direc tion of slip on the f ault p lane given the stress tensor follows from the 
formalism of iRamsev and Lisld (Il983h 

(8) l{m 2 (j, + n 2 ) : m{n 2 {\ -(f)- l 2 (f>) : -n{l 2 + m 2 (l - 0)) 

where I, m, n, are the direction cosines of the fault plane in the stress tensor 
reference frame and <j> is the stress shape ratio 

(9) 

ax — (T3 

This direction of slip, as the normal to the auxiliary plane, completes the definition 
of the focal mechanism. 

3.2. Generating error. I used a Fi sher distributio n to generate random directions 
and angles to perturb orientations ( Mardial . 1 19721 p. 231), and used quaternions 
to carry out the required rotations. Figure [T] shows an example of such Fisher 
distributed error angles and orientations plotted on a polar Schmidt net. Values of 
K for the Fisher distribution which correspond to perturbations of 1°, 5°, 10°, 15° 
and 20° can be obtained through interpolation from equation @] 

There are two ways in which to incorporate errors in the synthetic data set. 
There is the natural way of perturbing each focal mechanism directly. It is also 
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(a) k = 263 (5°) (b)K = 66(10°) (c)k=17(20°) 



Figure 1 . Polar Schmidt plots of examples of Fisher distributed 
errors, as used to perturb the generating stress tensor. The crosses 
are 10° apart 



possible to perturb the generating stress tensor before calculating the slip direction 
on each fault plane in the data set. The assumptions behind these two ways are very 
different, and I discuss their implications, statistical as well as physical, following 
the simulation results. 

3.3. Data inversion and statistics. Each synthetic data set is then processed 
by FMSI, which inverts the focal mechanisms to obtain the stress tensor. Thanks 
to FMSIETAB, part of the FMSI suite of programs, it is possible to generate a list 
of angular deviation between each fault plane and any given stress tensor. Such 
lists, comparing the deviations between focal mechanisms with the generating stress 
tensor as well as with the inverted stress tensor allow the calculation of the respec- 
tive spherical mean R and the concentration parameter k, using equations Q] and [5] 
Each combination of N fault planes with the different amounts of perturbation is 

re plicated 50 times. 

Hardebeck and Haukssonl ( 2001 ) proposed that for appropriate confidence in- 



tervals the correct stress tensor should fall within the P confidence region for an 
approximate proportion P of all the data sets. A plot of (sorted) probabilities ver- 
sus proportion should be a diagonal. If the confidence intervals are too wide, the 
plot will show a curve above the diagonal and too narrow confidence intervals will 
result in a curve below the diagonal. Thes e grap hs are an equivalent of P-P plots 



result m a curve below the diagonal. Ihe sc graphs are an equivalent ol Jr-Jr plots 
( Wilk and Gnanadesikanl Il968l : Holmgrenl Il995l ). P-P plots are scatter plots of 



Fi{qi),F 2 (qi), where F 1 (q i ) is obtained by applying 

(io) F 1 (F^ 1 ( Pi )) 

to the two empirical cumulative density functions (F) of the two data sets being 
compared. Here, I show P-P plots which compare the cumulative density function 
of the misfit against the inverted stress tensor with the cumulative density function 
of the misfi t agains t the generating stress tensor. 

Gephart ( 1990b ) proposed a modification of the method proposed by IZizica"sl 



(1955) to represent the orientation of a fault plane with respect to the principal 
stress directions. Gephart's unsealed Mohr Sphere uses the stress components a, 
ts and T6 as axes. 

The maximum shear stress and mean normal stress are respectively 

01 — 03 0i +03 

(H) T m = ^ , 0m = 7, 
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Figure 2. P-P plot of k for iV=20, with the generating tensor 
values along the x-axis and the inverted tensor values along the 
y-axis. The amount of perturbation in degrees is shown on each 
individual curve. 



The absolute values of the stress tensor are inaccessible, but the relative magnitudes 
of the principal stress components can be calculated. Their relation is given by 

(12) 0=^1 

03 - CT 1 

The three normalized stress components acting on a fault plane can be described 
as 

_/ _ _ _/ _/ 

ni\ ~ 11 m ^ 12 _ °i3 

and are completely determined by the dimensionless quantities (f> and /3y, where 
fiij is determined from 

(14) 4 = aufckPjl 

which relates the stress tensor components between the coordinate systems of the 
fault plane and of the regional stress tensor. The shea r stress will match the slip 



direction on a plane when a[ 2 is zero (|Gephartl . Il990b| ) 



Mohr Sphere plots of fault planes are often illuminating and are of great as- 
sistance with the interpretation of the inversions. This is also the case with the 
differently treated synthetic data sets, and I include a number of such Mohr Sphere 
plots for some of the data used. 

4. Results 

4.1. P-P Plots. The P-P plots obtained from the new analyses presented here 
(Figure [2HU) show a considerable improveme n t for the method compared to the 
results obtained by lHardebeck and Hauksson ( 2001 ). It appears that the change 



from linear statistics to directional statistics improves the appropriateness of the 
confidence interval estimation significantly. The graphs show a number of trends, 
including some unexpected ones. 

The most obvious, and expected, trend is a decrease in precision of the confidence 
interval estimates as the amount of error increases. This trend is present in each 
individual graph (Figures [2] [4| . 
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(a) Focal mechanism perturbation (b) Stress tensor perturbation 

Figure 3. P-P plot of k for iV=50, with the generating tensor 
values along the x-axis and the inverted tensor values along the 
y-axis. The amount of perturbation in degrees is shown on each 
individual curve. 




Figure 4. P-P plot of K for iV=100, with the generating tensor 
values along the x-axis and the inverted tensor values along the 
y-axis. The amount of perturbation in degrees is shown on each 
individual curve. 



What is unexpected is that if the estimate of the confidence intervals is im- 
precise, it is systematically too large. In the face of error, the method performs 
be tter than one would expect, som ething that was very pronounced in the study 
bv lHardebeck and Haukssor] ( 2001 ). 

The other noticeable and expected trend is an increase in precision of the confi- 
dence interval estimates as the number of focal mechanisms increases. This trend is 
very pronounced in graphs of the stress tensor perturbation set (Figures [5b, [3b and 
3]d). The graphs of the focal mechanism perturbation set (Figures [5^,, [3^ andQ]}) 
show an increase in precision for the smaller amounts of error, but an unexpected 
decrease for the larger amounts of error. 

There is also a clear tendency for an increased precision for the stress tensor 
perturbation sets against their focal mechanism perturbation equivalents (the (b) 
subfigure against the (a) subfigure for each of Figures [2j [3J and [4} . 
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(a) Generating Stress Tensor 




(b) Inverted Stress Tensor 



Figure 5. Mohr Sphere projections of the poles of the fault planes 
relative to the stress components defined by the respective stress 
tensor. The fault planes have been subjected to a 1° Fisher dis- 
tributed perturbation. 




(a) Generating Stress Tensor 




(b) Inverted Stress Tensor 



Figure 6. Mohr Sphere projections of the poles of the fault planes 
relative to the stress components defined by the respective stress 
tensor. The stress tensor has been subjected to a 1° Fisher dis- 
tributed perturbation. 

4.2. Mohr Sphere Projections. The Mohr Sphere projections (Figures [5H5]) arc 
examples of 4 data set of 50 focal mechanisms which have undergone cither a 1° 
(Figurcs[5]and[H]) or a 15° (Figures[7]and|H]) Fisher distributed amount of error. Each 



8 



S. A. REVETS 




(b) Inverted Stress Tensor 



Figure 7. Mohr Sphere projections of the poles of the fault planes 
relative to the stress components defined by the respective stress 
tensor. The fault planes have been subjected to a 15° Fisher dis- 
tributed perturbation. 



figure shows two sets of Mohr Sphere plots: one illustrating the fault planes against 
the generating stress tensor (or its spherical average when it was the generating 
stress tensor that underwent perturbation), and one showing the same set of fault 
planes against the inverted stress tensor calculated by FMSI. 

The first trend which stands out is the larger amount of scatter on the Mohr 
Sphere projections for the fault plane perturbed data sets (Figures [5] and [7]) com- 
pared to the projections for the stress tensor perturbed data sets (Figures |6] and 
HJ which is particularly noticeable on the a — Tb graphs. This increase in scatter 
includes the presence of fault planes with the wrong direction of slip (t$ < 0, seen 
as crosses which plot below the X-axis in the a — t$ and Tf, — rg graphs) . 

A second trend is the larger amount of change between the Mohr Sphere plots 
(going from the plot of the fault planes in the generating stress tensor Mohr Sphere 
to the Inverted tensor Mohr Sphere) for the data sets that have undergone fault 
plane perturbations compared to the data sets that have undergone stress tensor 
perturbation. There is more change visible going from the (a) subfigures to the (b) 
subfigures in Figures [5] and [7] than in Figures [6] and [8] 



5. Discussion 



Theoretical and statistical reasoning requires us to use directional statistics in- 
stead of linear statistics when we deal with focal mechanisms and deviatoric stress 
tensor inversion. The Fisher distribution and its properties are well established and 
its use in the present, theoretical study has been effective. The question does arise 
if this distribution is a good representation of actual seismological data. 

Discussions of the distribution of data misfit to stress tens o rs in the lit e rature 
(|Gephart and ForsvthL 1 1984 lHardebeck and Hauksso"nl l200ll IWvss et all Il992h 
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(a) Generating Stress Tensor 




(b) Inverted Stress Tensor 



Figure 8. Mohr Sphere projections of the poles of the fault planes 
relative to the stress components defined by the respective stress 
tensor. The stress tensor has been subjected to a 15° Fisher dis- 
tributed perturbation. 



consistently mention the non-normality of the data and refer qualitatively to expo- 
nential or Poisson distributions. Turning to spherical statistics resolves this issue 
immediately. The Fisher distribution, as gi ven in equation [51 shows directly its 
exponential nature. The method adopted by lHardebeck and Haukssonl ( 200ll) for 
their simulations, in which they combine exponentially distributed dip angles with 
uniformly distributed azimuth angles, generates in effect a set of Fisher distributed 
orientations. 

Since we are trying to assess the confidence intervals of the stress tensor inver- 
sion through simulations, we have to ensure that we use appropriate models and 
calculations to compare the perturbation distributions. It is highly instructive to 
compare the effects of injecting uncertainty into the simulations at two different 
points. It is for this reason that I opted to inject uncertainty or error in the fo- 
cal mechanisms to test the inversion process by perturbing the orientation of the 
generating stress tensor, as well as the more intuitive approach of perturbing the 
orientation of the individual focal mechanisms directly. 

The fault planes are chosen from uniform spherical random orientations (equa- 
tion [7J and the direction of slip is calculated from the stress tensor according to 
equation [8] It is clear that the relation between orientation of any of the focal 
mechanisms and the stress tensor is highly non-linear. That non-linearity neces- 
sarily also applies to any perturbation and in particular to the way in which it 
propagates in any calculation or inversion. It is incorrect to assume that a per- 
turbation distribution applied to a set of focal mechanisms would translate to the 
same perturbation distribution of the inverted stress tensor. This effect can be 
clarly seen in the Mohr Sphere plots where the scatter is much larger in the plots 
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of the data sets subjected to fault plane perturbation (Figures [5] and [7| than those 
subjected to tensor perturbation (Figures [6] and [8} • 

Comparing the dispersions of the data misfit to the inverted stress tensor with 
the dispersion of the perturbations in this case will be misleading. The dispersion 
of the stress tensor orientations (consistent with the perturbed focal mechanisms) 
is necessarily larger than the dispersion used perturb the focal mechanisms. It 
is therefore not surprising that the inversion method will converge to the correct 
stress tensor more often than one would expect from the dispersion used. This 
problem can be avoided by letting the perturbation occur on the generating tensor, 
before generating the focal mechanisms. This effect is clearly shown by the contrast 
between the P-P plots of the data sets subjected to fault plane perturbation and 
the data sets subjected to stress tensor perturbation (respectively the (a) and (b) 
subplots in Figures SHU • 

Perturbing the fault planes directly in effect introduces bias into the data set 
because each of the perturbed fault planes appears to have been generated by 
a potentially very different stress tensor, including incompatibly oriented stress 
tensors. This bias contaminates the dispersion. Using this contaminated dispersion 
to measure the accuracy of the confidence intervals then leads to the erroneous 
conclusion that the inversion method overestimates the confidence intervals. 

As far as simulations are concerned, the statistics can only legitimately be com- 
pared if they are commensurable. That is not the case when error is injected at 
the level of the individual fault plane. Unfortunately, this situation mirrors ex- 
actly what happens in the real world: the parameters describing a calculated focal 
mechanism are subject to error and uncertainty, not the generating stress tensor. 
Therefore, it may seem that the fact that FMSI yields correct confidence intervals 
only when the generating stress tensor is subjected to error, but fails to do so when 
focal mechanisms are error prone, is of no practical benefit. 

The paradox can be resolved as follows. I have demonstrated that FMSI yields 
correct confidence intervals when the dispersions used to calculate the confidence 
intervals are commensurable. I have also shown that error in focal mechanisms 
introduces bias in the statistics of the population at hand. But thanks to the 
side-by-side simulations, it is possible to estimate to what extent this bias contami- 
nates the calcu lations of FMSI. As one would expect from the law of large numbers 
(|ReveszL[l968h . a statistic converges to its theoretical value as the population size 
increases. This can be seen very clearly in the P-P plots: the results of the simula- 
tions converge to the diagonal as the sample size increases. This is very obvious for 
the simulations in which the generating stress tensor was subjected to error (the 
(b) subplots in Figures EHIJ . Closer scrutiny of the P-P plots of the simulations in 
which the focal mechanisms were subjected to perturbation shows a similar, albeit 
slower convergence (the (a) subplots in Figures EHIJ . There is a trade-off between 
amount of error and population size: FMSI is capable of estimating confidence 
intervals correctly up to a certain amount of bias. Figure [2^i shows that for pop- 
ulations of 20 focal mechanisms, FMSI is capable of dealing with an average of 5° 
of error on the focal mechanisms. This increases to 10° for 50 focal mechanisms 
(Figured) and to 12-13° for 100 focal mechanisms (Figured^,). 

The calculations by FMSI do not use directional statistics and the questions 
arises if a simple conversion is possible from the dispersion measure given by FMSI 
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Figure 9. Empirical relation between the R values and mean ab- 
solute deviation (°) for the Fisher Distribution 

to directional statistics and allow the application of the correct calculation of confi- 
dence intervals as shown here. The mean deviation calculated by FMSI is an average 
of angles while the circular mean is the average of the cosines of angles. There is 
no analytical formula to go from one to the other. Nevertheless, there is a clear 
relation between the mean absolute deviation and R for the Fisher Distribution 
(see Figure IHl) and graphs of this nature could be used to estimate values which can 
then be used to determine the confidence intervals. The correct procedure, albeit 
more calculation intensive, is to use the list of individual focal mechanism misfits 
(which can be obtained through FMSIETAB) and to apply equation [S] 

6. Conclusions 

The processing and discussion of data from focal mechanisms and deviatoric 
stress tensors should use directional statistics. 

The genesis of synthetic data in the context of stress tensor inversion requires 
very careful scrutiny of how perturbation is incorporated. 

The inversion of deviatoric stress tensors from focal mechanisms using the grid 
search method of Gephart yields reliable confidence intervals if correct, directional 
statistics are used. When calculated misfits are small for moderate to large numbers 
of focal mechanisms, the method recovers the true stress tensor and confidence 
intervals become superfluous. 

Gephart 's FMSI is a highly effective and reliable method for stress tensor inver- 
sion. 



7. Data and Resources 

The synthetic data, data inversion and additional calculations relied on a com- 
bination of_B ash shell scripts and procedures written for the Octave program 
( Eatonl 2002). The FMSI suite of ] urograms is in the public domain and made 
1990a|) . The fortran source code is freely available from 
www.geo.cornell.edu/pub/FMSI. 
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